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We propose a new concept for the design of high-efficiency photocells based on ultra-thin (sub- 
micron) semiconductor films of controlled thickness. Using a microscopic model of a thin dielectric 
layer interacting with incident electromagnetic radiation we evaluate the efficiency of conversion 
of solar radiation into the electric power. We determine the optimal range of parameters which 
maximize the efficiency of such photovoltaic element. 



Improvement of efficiency of semiconductor photo- 
voltaic elements (solar cells) has been an important tech- 
nological challenge for several decades. The maximum 
possible efficiency obtains when every incident photon 
generates an electron- hole pair, which then separates into 
electron flowing to cathode and hole flowing to anode 
The limitations that reduce the efficiency of the practi- 
cal solar cells relative to the ideal are 1) light reflection 
at the interfaces, 2) incomplete absorption of light en- 
tering the device due to finite thickness, 3) electron-hole 
relaxation inside the absorbing medium during diffusion 
to the leads [3, 0). Interplay between two latter mecha- 
nisms leads to an existence of an optimal device thick- 
ness, typically a few optical wavelengths. Here we show 
that the interface reflection, commonly considered a com- 
pletely independent loss mechanism, shows an interest- 
ing interplay with absorption in ultra-thin film devices. 
This opens a possibility for a new generation of ultra-thin 
(sub-wavelength) photovoltaic elements with efficiencies 
rivaling the best conventional devices. 

A "working body" of a solar cell is typically a semicon- 
ductor with relatively high absorption index at frequen- 
cies corresponding to those of the sun quanta fko sun ~ 
fesTsunj Tsun — 6000 K. Such semiconductors, however, 
process a rather high refraction index n at these frequen- 
cies. As a consequence, a fraction (n — l) 2 /(n + l) 2 of 
the incident light is reflected from the surface of the think 
device. To reduce this loss, often anti-reflective coating 
are applied to the surface of the device. On the other 
hand, for sub-wavelength thin films, the reflection can 
be significantly smaller (for a reason similar to why even 
metallic films are transparent when thin enough). Thus 
reducing the film thickness one should reach an optimum 
where reflection is reduced but the absorption is still sig- 
nificant. 

Also, in such thin devices the carrier recombination 
is naturally reduced. Electron-hole recombination which 
prevents efficient charge separation in the photocell is a 
major limiting factor in device operation. There are nu- 
merous mechanisms which lead to the charge relaxation 
in a bulk of a semiconductor. These mechanisms include 
spontaneous emission as well as phonon or impurity in- 
duced relaxation. While it is difficult to control these 
processes in the bulk of a semiconductor, it is clear that 



their contribution can be significantly reduced if diffu- 
sion length of electrons and holes is large compared to 
the width of the semiconducting layer. 
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FIG. 1: Insets (a) and (b): Schematics of the device. Insets 
(c) and (d): Band structure of the device without and with 
external load. 

Diffusion length for most semiconductors used in pho- 
tocells is of the order of a few microns. Since this distance 
is comparable to a typical wavelength of the sunlight, 
one could expect that the specific absorption (the ratio 
of the absorbed power to the incident power of radiation) 
in such a thin semiconducting layer is insufficient for any 
practical use. In order to see whether this is the case, 
it is instructive to look at the absorption of radiation 
in a layer of thickness d, e.g. Fig. l(a,b). For simplicity 
we assume that the radiation is incident perpendicular to 
the surface of the layer and is monochromatic with wave- 
length A. The layer has a dielectric constant whose real 
and imaginary parts are CR e and ei m respectively. The 
specific absorption of the dielectric layer can be easily 
evaluated by solving the wave equation (n/c 2 )A = d 2 A 
(A is radiation field vector-potential, n — e 1 / 2 is the re- 
fraction index, and c is the speed of light) in the three 
regions, i.g., Fig. 1(b), and taking into account the conti- 
nuity conditions at the boundaries of the dielectric layer, 
Ai = A2, A2 — A3, d z A\ = d z A2, etc. After straight- 
forward algebra one finds that the specific absorption is 
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^abs = 1 — \ t\ 2 — \ r \ 2 , where the amplitudes of transmitted 
and reflected waves are 

, 4n exp (id/X) 

(n + l) 2 exp (-ind/X) - (n - l) 2 exp (ind/A) ' 1 &) 

(n 2 — l)[exp(md/A) — exp (— ind/X)] . 

r = -(lb) 

(n + l) 2 exp (—ind/X) — (n — f ) 2 exp (ind/X) 

The specific absorption evaluated according to the above 
equations is presented in Fig. 2 for a GaAs slab as a 
function of its thickness d. GaAs has a relatively nar- 
row (~ 1.4 eV) bandgap and therefore is widely used 
in high-efficiency photocells. Since dielectric function 
of GaAs is strongly frequency dependent [4], in Fig. 2 
we plot the specific absorption for several energies typ- 
ical to the quanta of solar radiation. The 2 eV curve 
corresponds to the relatively low imaginary part of the 
dielectric constant and thus saturates slowly exhibiting 
several oscillations due to the interference between re- 
flected and transmitted components. The 3 eV and 4 eV 
curves correspond to much higher absorption (for exam- 
ple e^ As (3 eV) ~ 17) and saturate much faster. Prior 
to saturation both curves exhibit a peak (again due to 
the interference) at roughly d ~ A/|e|. Remarkably the 
value of the specific absorption at the peak (~ 0.42 at 
d ~ \/\e\ ~ 20 nm for 3 eV curve) is nearly the same as 
its saturation value (0.51 at d ^> A). Thus we conclude 
that the solar radiation can be absorbed by a semicon- 
ducting layer of submicron thickness almost as efficiently 
as by an infinitely thick slab. 

In this paper, following the simple above considera- 
tions, we propose a new concept for the design of pho- 
tovoltaic elements based on thin semiconductor films of 
controlled thickness. To put our arguments on more rig- 
orous footings in the following we consider a detailed 
microscopic model of a dielectric layer interacting with 
solar radiation. Effective Hamiltonian can be derived 
from standard quantum-mechanical interaction between 
matter and radiation, (e/mc)A(r)-p + (e 2 /2mc 2 )A 2 (r), 
where here and in the following we assume Coulomb 
gauge for the electromagnetic field. The radiation in- 
duces transitions between valence and conduction band 
of the semiconductor. We assume that the temperature 
of the semiconductor is 0, and therefore these are the 
only possible transitions in the system (the valence band 
is completely full and the conduction band is empty) . De- 
noting the Bloch states for the valence(conduction) bands 

as Vk (c) (r) - ex P HkrK (c) (r), = E^^ 

(Ho is a Hamiltonian of the crystal in the absence of 
coupling to radiation field) one can rewrite the radiation- 
matter interaction Hamiltonian in terms of single particle 
states in the semiconductor as 

H ^ = ^c E Kkl«S a >44- qia C+H.c.. (2) 

We make the following assumptions: (1) While the 
electromagnetic field does not significantly vary with 



distance inside the film, the electronic wave-functions 
are effectively 3-dimensional - we assume that A sun ~ 
2irhc/(k B T sun ) > d > h/ (m* Eg) 1 / 2 , where d is the 
thickness of the film, to* is exciton effective mass, 
(to*) -1 = to" 1 + to^ 1 , to„( c ) are effective masses in va- 
lence and conduction bands, and E g is the band-gap (in 
this paper we assume that bands have extremuma at zero 
momentum). This assumption allows one to carry out an 
analytic calculation with rather simple and transparent 
results. We will discuss the validity of this approxima- 
tion at the end of the paper; (2) The bands have different 
symmetry, say s and p, so index a denotes angular mo- 
mentum of an electron in p band. Since wave-vectors 
of the incident radiation are nearly perpendicular to the 
surface of the film and the film is assumed infinite in 
x — y dimension, only in-plane components of the an- 
gular momentum (a = x,y) are relevant in Eq. ([2]); (3) 
The coupling in Eq. ((2|) is isotropic and the coupling con- 
stant t a = (uo\ep a / (mc)\uo a ) ~ Egpj '(cS 1 / 2 ), where p is 
the effective dipole moment per unit cell of the film and 
S is the surface area of the film. (4) Due to external 
electric load the bands have effectively different chem- 
ical potentials, p n and fx p , e.g., Figs. 1(c) and 1(d). 
That is, once electron is promoted from valence to con- 
duction band, it immediately "rolls over" to the left lead, 
which corresponds to the infinite transition rate between 
the semiconductor and the metallic lead (an infinitely 
thin Shotkey barrier) . Clearly, were the rate comparable 
or slower than the electron-hole relaxation rate, the effi- 
ciency of the cell would have decreased. (5) The p — n 
junction is prepared (doped) as shown in Fig. 1(c) - the 
top of the valence band in the p-doped area of the junc- 
tion lies just below the bottom of the of the conduction 
band in the n-doped area. Therefore the maximum volt- 
age the photo-cell can sustain is equal to E g /e, which 
corresponds to the assumption of maximum efficiency of 
Shockley and Quasser[lj, i.e., each electronic transition 
from valence to conduction band generates energy E g in 
the circuit. The photocurrent is defined as the rate of 
the charge transfer between the valence and conduction 
band, I p h = [H m t , Ylk c ik Ck ] ■ To the lowest non- vanishing 
order it can be expressed as 

J P h = J £ ~ K-^z = z' = 0). (3) 

k qj_a 

The photocurrent of Eq. ([3} is independent of the 
voltage across the cell as far as it does not exceed 
the bandgap of the semiconductor, e.g., Fig. 1(d). 
For larger bias, within our assumptions, a reverse 
current begins to flow. In Eq. ([3]) D<„ (a;, z, z') = 
J dtd 2 r ± expi(ujt + r ± ■ q x ) (T K A~ (t, r ± , z)A+(0, 0, z')) 
is the "lesser" Green's function of the electromagnetic 
field defined along the Keldysh contour, where super- 
scripts ± denote forward and return branches of the 
contour The Green's function is inhomogeneous 
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FIG. 2: Specific absorption of GaAs slab as a function of its 
thickness. Different curves correspond to different frequencies 
(energies) of incident radiation. 



along z-direction, i.e., perpendicular to the film surface. 

In order to incorporate effects of absorption and reflec- 
tion from the film, it is necessary to include renormaliza- 
tion of the photon Green's function due to interactions 
with the film according to Eq. ([2]). These effects can be 
treated by means of Dyson equation, which, for the case 
of two-dimensional film reads 

D a /3{uJ,q±,z,z') = Doa/3(w,q_i_,z- z') 
+Ooa 7 (w, q±, z)t, 7 s(u, qx)D S 0(v, q_L, 0, z') , (4) 

where hats denote the standard 2x2 matrix struc- 
ture of the non-equilibrium Green's functions. The 
self-energy E 7d - in Eq. (J4J is quasi- two-dimensional. 
Since I h/(m*E g ) 1 / 2 , effects related to the finite 
width of the semiconducting slab can be neglected and 
S 75 (w,qj_,0) = J2 q , £>~fs(u,<l), where S 7(5 (w,q) is the 
self-energy defined for the bulk of the semiconductor. 
Moreover, due to x — y symmetry S 7 5 reduces to a 5 7 aS, 
where £ depends only on |q_i_|. Therefore we obtain a 
closed form equation for D a/ 3q ± (u>, z = z' = 0): 

D a ^(0,0) = J D 0a/3 (0) + J D 0a7 (0)X;£) 7 ^(0,0), (5) 

where components of the bare Green's function -Do a/3 f°r 
solar radiation are: 



D, 



R(A) 
Oq/3 



Kqx,0) 



^ 47r (^a/3 ~ QaQ/3/q ) 
1* 



■ wl ± iS 



X [S(iU — LUq)flq — 8(b) + Wq)(l + U-q)] . 



(6a) 



(6b) 



In Eq. 



hc\q\ and the n q is the distribution 



nq = n*9(q z )6(4>q z - |qj_|), where is Bose 



radiation wavevectors are uniformly distributed within a 
cone with an opening angle 2(f) (4> <C 1). Moreover, in 
order for the incident power to be maximum, we assume 
that the surface of the cell is perpendicular to the cone's 
axis. Then 

distribution function with temperature T E 

Furthermore due to the homogeneity of the Green's 
function D a pq ± in x — y plane one can seek for solution 
ofEq. §5§ in the form D a p — £>i5 a p + D 2 q a qp/q 2 , where 
£>H2) are 2x2 matrices in the Keldysh space, but depend 
only on the absolute value of the wavevector q. Substi- 
tution of this ansatz into Eq. ([5]) yields two independent 
equations for D\ and for D3 = D\ + £>2- After solving 
those equations one finds 
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01(3) J 



1 -^01(3)('^C)1(3)/ 



(3); 



E«][(L> ( A 



01(3)' 



(7) 



where £>oi is the diagonal part of Doap, e.g., Eqs. ([6]), 
and D03 = -D01 + ^02, where -D02 is the transverse part 
of Doafj. E fl ( A ' and £< are retarded(advanced) and 
"lesser" parts of photon self-energy. Also we find a stan- 
dard expression for the retarded(advanced) Green's func- 
tions of the radiation: 
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R(A) 
1(3) 



LV- L/ 01(3)' 1 



(8) 



function of solar radiation. We assume that the incident 



We can now evaluate the photo-current in Eq. ([3]) in 
terms of the Green's function of Eq. . The first contri- 
bution is due to absorption of incident radiation accom- 
panied by transfer of electrons from valence to conduction 
band. It comes from the first term in the numerator in 
the RHS of Eq. |J7J| . This term can also lead to the reverse 
current due to spontaneous and stimulated emission, i.e., 
transitions from conduction band to the valence band ac- 
companied by creation of a real photon. It arises due to 
the S(uj + uiq) term in Eq. (|6b[) . This process is, however, 
not allowed while the energy gap E g exceeds the applied 
voltage jtt„ — fj,p, e.g., Fig. l(c,d). In this situation the 
cell becomes a light-emitting diode, and, as was stated 
above, we are not interested in such case in this paper. 
The second contribution comes from the £ < term in the 
RHS of Eq. {7}. It corresponds to the emission of vir- 
tual quanta of radiation by one electron-hole pair and 
their subsequent re-absorbtion by another pair, resulting 
in an incoherent simultaneous transfer of two electrons 
from conduction to valence band. This process gives a 
reverse contribution to the current which, once again, is 
non-zero only in the "light-emitting" regime (or at high 
temperature). 

The self-energies in Eqs. (|7|8|) can be evaluated in 
terms of the electronic Green's functions. The leading 
contribution comes from the conventional polarization di- 
agram, i.e., a convolution of two electronic Green's func- 
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tions. For non-equilibrium situation one obtains 



+ G&(wW)G&(w , ) + ('>~<0 



(9) 



and E A = In Eq. © Gf ( ( c j } are retarded (ad- 

vanced) Green's functions of valence(conduction) elec- 
trons and G^ c ^ is the Keldysh Green's function. Note 
that in Eq. © we evaluated the self-energy at zero 
wavevector, since photon wavevectors are small com- 
pared to those of electrons, and therefore the self-energy 
£ is weakly dependent on q for direct bandgap materials. 

The self-energies can be easily evaluated for non- 

interacting electrons. ThenG^f^(w) = {lu— E^^ztiS)^ 1 

G* c)k H = (1 - 2< (p £)J( W - E<\ where < (p)k 
are Fermi filling factors. The valence and conduction 
electrons are assumed to have chemical potentials corre- 
sponding to those of the two leads, y, n and fi p respec- 
tively. For these Green's functions the imaginary part of 
the self-energy yields: 



aE, 



1/2 



J^ A \ Ut 0)=±^-Q(u;-E B )y/^E^, (10) 

Z7T 

where we have introduced a dimcnsionless parameter 
a = p 2 dEg /2 (m*)^ 2 /{2 1 / 2 h i c). A similar calculation 
for the real part of yields an ultraviolet diver- 

gence. This, however, is an artifact of our approxima- 
tion of the infinitely thin absorbing layer (similar prob- 
lem occurs in the quantum electrodynamics treatment 
of electron-photon interaction) . In a proper microscopic 
theory this divergence in the real part of the self-energy 
is exactly cancelled by the e 2 A 2 (r)/(2mc 2 ) term in the 
radiation-matter interaction Hamiltonian (the frequency 
sum rule). Taking into account this cancellation is equiv- 
alent to performing the Kramers-Kronig transformation 
on the dielectric function [Re e(ui) = 4-7rc 2 Re £(w)/oj 2 ], 
rather than self-energy: 



7T J UJ' (Id' 



(u> - to) 



aE. 



1/2 



2tt 



E g \~ 



~E g \). 



(11) 



Note that in Eqs. (|10|lip we used the three-dimensional 
expression for the self-energy and therefore we can re- 
express parameter a in terms of the conventional zero 
frequency dielectric constant eo, a = 2lE g (eQ — l)/(hc). 

From Eqs. (|6I7I10I11|) one can evaluate the photocur- 
rent 7 p h given by Eq. ([3]). In this paper we are interested 
in the maximum efficiency of the photocell, which can 
be defined as r? max = (P7) max /P in , where (IV) max is the 
power that is dissipated in the circuit assuming an opti- 
mal load and P; n is the power of incident solar radiation. 




FIG. 3: Dependence of photovoltaic efficiency on dimension- 
less parameters a and 6. 



Since the photo-current is only weakly dependent on volt- 
age for V < Eg and becomes negative due to spontaneous 
and stimulated emission emission for V > E g , we have 
(7V) max ~ IphEg/e. The incident power per unit area 
P; n = cu, where energy density u = (2/v) X^q^q^q! 
where v is the mode quantization volume and factor 2 is 
due to two polarizations of the light waves. Carrying out 
the calculation we obtain the following closed form ex- 
pression for the maximum photovoltaic efficiency of the 
cell: 



_ 30ab 4 

7T 



dx 



cxp (bx) — 1 



x' J \fx 



,(12) 



(x + aVx - l) 2 + a 2 (2 - Vx + 1 - \J x - 1) ; 

where we introduced another dimensionless parameter 
b = E g /T sun . Function rj max (a,b) shown in Figure 3 has 
a pronounced maximum at a ~ 2.0 and b ~ 2.4. This 
maximum corresponds to the first maximum of the spe- 
cific absorption curves in Fig. 2. Note that since our 
microscopic theory is valid only for thin dialectic layers 
(d <C A sun ), it does not account for the saturation of 
specific absorption at larger thickness, i.e., when d > A. 
However, since r) mayi reaches maximum at d ~ A sun /eo, 
our theory is fully self-consistent for semiconductors with 
sufficiently high value of the dielectric constant. 

The optimal value of parameter b corresponds to the 
bandgap energy E° pt ~ 1.2 eV. Since this value is close 
to the bandgap energy in GaAs (P^ aAs ~ lAeV), we 
conclude that GaAs is a good candidate for the prac- 
tical realization of thin film photocells. According to 
Eq. (fT2|) the optimal thickness of such GaAs layer is 
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d opt ~ 1.2tic/[E°P t (e^ aAs - 1)] ~ 15 nm, which is within 
the fabrication capabilities of contemporary molecular 
beam epitaxy technology. Another promising material 
is amorphous Silicone which, unlike crystalline Si, has 
large imaginary part of dielectric constant Q. Moreover, 
in ultra-thin devices considered here, the thermal equi- 
libration of photo-generated carriers may occur on the 
time scales longer than the charge separation timescale. 
Thus the hot-carrier physics may lead to further en- 
hancement of the efficiency. 
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